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Abstract: As the frontline encountering various microbes, intestine plays critical roles in insect immune 
responses. In order to investigate the immune responses in the silkworm midgut, we identified 18 
differentially expressed genes after oral infection by Pseudomonas aeruginosa and Staphylococcus aureus 
using annealing control primer ( ACP) based reverse transcriptional PCR technique. The expression 
profile of four genes during the first 24 hours post infection was monitored by quantitative real-time PCR. 
The results showed that the peptidoglycan recognition protein L1 (PGRP-L1) gene and a serine protease 
precursor gene were up-regulated by P. aeruginosa infection exclusively, while the 30kP protease A 
precursor gene was up-regulated significantly by both P. aeruginosa and S. aureus infection. Our study 
identified genes involved in recognition of invading bacteria and immune signaling pathways in the larval 
midgut of the silkworm after oral infection and may give some clues for further functional investigation of 
these genes. 
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1 INTRODUCTION 


Insects have evolved elaborative and effective 
immune systems to deal with encountered microbes. 
By introducing bacteria into Drosophila body cavity 
with injection or pricking, it has been established 
that the production of antimicrobial peptides 
(AMPs) through Toll and IMD pathways functions 
as systemic immune reactions ( Ferrandon et al., 
2007; Vallet-Gely et al., 2008). Realizing ingestion 
as the main route for bacteria and other pathogens to 
enter into the hosts, researchers started to study the 
local immune defenses in the intestine after infection 
through the natural way, oral ingestion ( Basset et 
al., 2000; Schmidt et al., 2008). The generation of 
reactive oxygen species (ROS) is an immediate 
defensive reaction to ingested pathogens. ROS in 
Drosophila reached the highest level at 3 h after 
natural infection with bacterium Erwinia carotovora 
carotovora 15 (FEcc15) (Ha et al., 2005). ROS and 
nitric oxide are critical determinants of Plasmodium 
survival in Anopheles gambiae ( Marois, 2011). In 
order to combat ROS-resistant pathogens, their hosts 
produce AMPs via IMD pathway in the gut as the 


second line of defense (Ryu et al., 2006; Ryu et al., 
2010). 

Little is known about gut immune reactions in 
insects other than the fruit fly and mosquito. 
Silkworm, Bombyx mori, has received much 
attention as a lepidopteran model insect in recent 
years. The transcript levels of key genes of the major 
immune pathways in silkworm were profiled 
microarray after oral 
infection by bacterium Bacillus 


bombyseptieus ( Huang et al., 2009). Toll-9 


transcription was up-regulated in the silkworm gut 


systematically based on 
pathogenic 


after challenge with Escherichia coli or Beauveria 
( Wu et al., 2010). 


Transcription of serine protease inhibitor ( serpin ) 


bassiana in the diet 


was also up-regulated in the midgut when the 
silkworm was infected by B. mori cypovirus ( Wu et 
al., 2011). 

In this study, we used annealing control primer 
( ACP ) PCR strategy to 
expressed genes after the 
challenged with ingested bacteria and detected 


screen differentially 
silkworm had been 


several inducible genes using quantitative real-time 
PCR. The efficiency of ACP technique and possible 


functions of these genes identified in this study were 
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also discussed. 


2 MATERIALS AND METHODS 


2.1 Silkworm rearing and bacterial infection 
B. mori ( Nistari ) larvae were reared on 
artificial diet (Fukui and Co., Ltd. ) at 26°C with 
12 h light/12 h darkness photoperiod. When B. 
mori was at day 1 of the 5th instar, tetracycline (50 
wg/larva) was added to the diet block (5 mm’) to 
kill indigenous bacteria in the silkworm ( Hamamoto 
et al., 2004). Pseudomonas aeruginosa ( PAO1 ) 
and Staphylococcus aureus were cultured in Luria- 
Bertani medium at 37°C overnight. The cells were 
then collected by centrifugation at 5 000 g for 15 
min, washed with PBS three times and resuspended 
with PBS to Ag = 20. Twenty four hours after the 
50 pL of the 


preparations or PBS (as the control) were added in 


antibiotic treatment, bacterial 
the diet block or injected into the hemocoel to 
challenge the larvae. Each group includes 20 larvae. 
The animal survival was recorded in the following 
five days. For RNA extraction, 24 hours after the 
insects finished the diet blocks, 10 larvae from each 
group were dissected and the midgut tissues were 
collected and stored in Trizol reagent ( Invitrogene ) 
at —80°C. 
2.2 Total RNA extraction and cDNA preparation 
The midgut tissues were homogenized using a 
pellet pestle ( Kontes) and total RNA was isolated 
using Trizol Reagent (Invitrogen) according to the 
standard protocol. The obtained RNA was purified 
using RNeasy MinElute Cleanup Kit ( Qiagen) and 
the remaining genomic DNA contamination was 
removed by DNase treatment ( Promega). Agarose 
gel (1. 0%) was used to check the total RNA 
quality by loading 1 uL of total RNA samples. The 
concentration and A,,,/A,., were determined by 


The purified RNA 


samples (3 wg) were used as the templates in the 


Biophotometer ( Eppendorf ). 


reverse transcription system of DEG GeneF ishing Kit 
(Seegene) to synthesize the first-strand cDNA using 
dT-ACP1 (1. 0 pmol/L ). 


transcription was carried out for 1 h at 42°C following 


primer Reverse 
the manufacturer’ s protocol. The first-strand cDNA 
was 5-fold diluted with water and stored at — 20°C 
for GeneF ishing PCR. 
2.3 ACP-based differential display RT-PCR 
Differentially expressed genes were screened by 
ACP-based PCR using the GeneFishing DEG kits 
according to the manufacturer’ s protocol. Standard 
two-stage PCR was carried out in a final reaction 
volume of 20 uL. In the first stage, the second- 


strand cDNA was synthesized using the following 
thermal cycling parameters: one cycle at 94°C for 5 
min, 50°C for 3 min and 72°C for 1 min. The 
second-stage PCR was carried out to amplify cDNA 
using the following program: 40 cycles of 94°C for 
40 s, 65°C for 40 s, and 72°C for 40 s, followed by 
extension at 72°C for 5 min. In our analysis, a total 
of 40 different ACP (annealing control primer) pairs 
were used. The amplified PCR products (5 uL) 
were separated on 2% agarose gels and stained with 
ethidium bromide. The remaining PCR products 
were stored at — 20°C for TA cloning. 
2.4 Cloning and sequencing of differentially 
expressed genes 

Bands with different 
between control and infected groups were recorded 


significantly intensity 
and the rest of the PCR products were separated on 
1. 5% low-melting-point agarose gels ( BioBasic 
Inc. ). QIAquick Gel Extraction Kit ( QIAGEN ) 
was used to extract PCR products from the gels. The 
recovered PCR fragments were ligated with the TA 
cloning vector pMD-19 (TaKaRa) and transformed 
into E. coli JM109. The transformed bacteria were 
cultured for plasmid extraction using Plasmid Mini 
Kit ( Omega ) and sequencing ( Shanghai Sunny 
Biotechnology Co. ). All sequences were subjected 
to BLAST search for gene identification at Silkworm 
Genome Database ( SilkDB: http://www. silkdb. 
org/silksoft/blast2-simple. html ), GenBank 
(NCBI) and KAIKObase (http://sgp. dna. affre. 
go. jp/KAIKObase ) . 
2.5 Quantitative real-time PCR confirmation of 
the selected genes 

To confirm the results of the ACP technology , 
quantitative real-time PCR was used to detect the 
amplification of four differentially expressed genes for 
A11, 30kP protease A precursor, PGRP-L1 and 
serine protease precursor. Total RNA was extracted 
from midguts of ten infected 5th instar larvae at 1, 
2, 8, 16 and 24 h post infection (hpi) as well as 
from the uninfected larvae by Trizol reagent 
(Invitrogen). The RNA was treated with DNase I 
and purified using Direct-zol™ RNA Miniprep 
(Zymo Research). M-MLV Reverse Transcriptase 
(Promega) was used in the reverse transcription 
reaction following the manufacturer’ s instructions. 
Real-time qPCR was performed using the SYBR 
Premix Ex Taq™ II Kit (TaKaRa) and Bio-Rad IQ5 
Real-Time PCR 
manufacturer’ s protocol. Diluted first-strand cDNA 
(2 uL, 50 ng/L) was used in each 25 uL reaction 


mixture and the final concentration of primers was 


system according to the 


0.4 pmol/L. A reference gene /F4A was used to 


5H 


normalize the samples ( Wu et al., 2010). Primers 
used in the real-time gPCR are displayed in Table 1. 
All the reactions were run in triplicate using the 
following thermal cycling parameters; denaturation at 
95°C. for 30 s, followed by 40 cycles of 95°C for 10 s 
and 53°C for 30 s. Following amplification, melting 
curves were constructed to ensure the veracity of the 
experimental results ( temperature range: 53 - 
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95°C). Standard curve of every gene was made to 
measure PCR efficiency, which was in the range 
from 0.9 to 1.1. R? values of the standard curves 
were greater than 0. 997. Relative quantitative 
method (AACt) was used to evaluate the expression 
difference. In each time point, the expression level 
relative to the expression level of the control, which 


was arbitrarily set to one, was shown. 


Table 1 Primers used in Quantitative real-time PCR for confirmation of differentially expressed genes 


Protein GenBank accession no. 
All BGIBMGA004115 
30 kP protease A precursor BGIBMGA012788 
PGRP-L1 AK378516 
Serine protease precursor BGIBMGA014427 
IF4A BGIBMGA003186 
3 RESULTS 


3.1 Mortality of silkworm to bacterial feeding 
and injection 

First we compared the silkworm mortality to 
different routes of infection. Direct injection of the 


a m z a a a 


Primer sequences (5’ —3’) Product size (bp) 


: TTCCTTACCTGGAGAATG 77 
: AACACAAATTGAAGTCACT 
: CTCTTACCATAAGGATTAC 161 
: CGATATTGTTGATTACGA 
: GTGAAGTTATCTACGGAATGG 152 
: ACAAATTGGGCTGACCTA 
: GGATTCAGAATTAACACTAACAA 161 
: GTCCGTGGATAATCATAGC 
: AATCAACGTGAGCGTGAAGTG 192 
: TCCCTTACGACCAAAACGTC 
bacteria into the hemocoel caused the death of most 
silkworms in the test groups in two days, with PAOI 
killing the worm faster than S. aureus (Fig.1:; A). 
Feeding infection showed slight pathogenic effect on 
the test silkworms (Fig. 1: B). 
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Survival rate of silkworm infected by injection (A) and feeding (B) infection 


Cont: Control (PBS); Pa; Pseudomonas aeruginosa; Sa; Staphylococcus aureus. 
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Fig. 1 
3.2 Identification of differentially expressed 


genes after bacterial infection 

To identify genes with a predominant change of 
expression in the silkworm midgut after natural 
infection, we examined the gene expression using 
GeneFishing, a PCR-based differential display 
technique. RNA from midgut tissues infected with 
bacteria was individually subjected to ACP-based 
RT-PCR using a combination of two anchor Oligo 


(dT) primers (dT-ACP1 and dT-ACP2) and 40 
arbitrary primers. PCR products were separated and 
visualized on 2% agarose gels. Several bands were 
defined as differentially expressed when there was a 
clear difference in intensity between infected and 
uninfected samples. A representative gel image is 
in Fig. 2. Totally, identified 18 
differentially expressed bands. 


shown we 
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Fig. 2 Electrophoretogram of PCR products from the differential gene expression analysis 
A2 — A109 are the primer pairs used in the GeneFishing PCR. Cont; Control (PBS) ; Pa; Pseudomonas aeruginosa infection; Sa; Staphylococcus aureus 


infection. The same for Fig. 3. 


The differentially expressed bands were cut and 
subjected to sequencing. The DNA sequence of each 
differentially expressed gene was analyzed by BLAST 
search in the SilkDB, GenBank and KAIKObase. 
Cloning, sequencing, and BLAST search revealed 21 
differentially genes (Table 2). Among them, A102c 
had no homolog in the databases; A11, A12 and 
A107b encoded proteins with unknown functions. Four 
genes encoding BGIBMGA004115, peptidoglycan 
recognition protein-L1 (PGRP-L1, AK378516) , 30kP 
protease A precursor ( BGIBMGA012788 ) and serine 
protease precursor ( BGIBMGA014427 ) , respectively, 
are of particular interest to us. 

3.3 Confirmation of the ACP observations by 
real-time qPCR 

Differential expression of the four selected genes 


was confirmed (Fig. 3). The expression of PGRP- 


Ll and serine protease genes was up-regulated 


exclusively by PAO1 infection. The expression of 
30kP protease A precursor encoding gene 
significantly increased 2 — 16 h post PAOI or S. 
aureus infection. BGIBMGA004115 responded to the 
infection quickly and dramatically. Its mRNA level 
increased about 5-fold at 1 h post infection, and then 
dropped about 1 000-fold during 2 —8 hpi compared 
to the non-infected samples. 


4 DISCUSSION 


The intestinal barrier plays critical roles in host 
immune system. When the silkworms were infected 
by gut or hemocoel insults, we saw different survival 
results. Hemocoel infection caused the silkworm to 
die quickly. This implies that most bacteria acquired 
from diet have been killed in the gut before they pass 
the midgut barrier and enter the hemocoel. The other 


5 H 


Table 2 Annotation of the differentially expressed genes by the Blast search 
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Clone SILK DB/GenBank/KAIKObase ID E-value Description of gene product cDNA 
A2 BGIBMGA004005 1.00E -108 Putative cuticle protein ( Japanese oak silkworm) AK378233 
A4 BGIBMGA012788 0 30kP protease A (43k peptide) precursor AK378510 
ASa BGIBMGA009130 1.00E -156 Ribosomal protein $20 AK386501 
A5b BGIBMGA004005 0 Putative cuticle protein ( Japanese oak silkworm) AK378233 
Aba BGIBMGA010979 0 Haemolymph juvenile hormone binding protein (JHBP) AK388525 
A6b BGIBMGA012414 3.00E -61 Ribosomal protein L17 AK384926 
All BGIBMGA004115 0 Unknown AK388527 
A12 BGIBMGA004115 1.00E -164 Unknown AK388527 
Al6a AB455925. 1 0 Transposable element; non-LTR_TREST-W AK381770 
Al6b NM_001099806. 1 1.00E -153 Small nuclear ribonucleoprotein polypeptide F AK385097 
Al6e AK378516 5. 00E -95 Peptidoglycan recognition protein L1 AK378516 
A17 AB473763. 1 0 Transposon mariner-like element, non-LTR retrotransposon LIBm AK378418 
A19 BGIBMGA009411 1. 00E -29 Ribosomal protein $16 AK383780 
Al02a BGIBMGA013756 0 Peritrophic membrane chitin binding protein AK378243 
A102b BGIBMGA002635 0 Sugar transporter AK383428 
A102c No hit - - No hit 
A104 BGIBMGA013413 1.00E -151 Tubulin alpha-1 chain AK388384 
A105 BGIBMGA010584 0 Chymotrypsinogen AK378306 
A106a BGIBMGA010979 1.00E -142 Haemolymph juvenile hormone binding protein ( JHBP) AK388525 
A106b BGIBMGA009477 1.00E -174  Carboxypeptidase A AK378467 
A106¢ BGIBMGA010979 0 Haemolymph juvenile hormone binding protein (JHBP) AK388525 
Al07a BGIBMGA009698 1.00E -125 Nascent polypeptide associated complex protein alpha subunit AK383828 
A107b AK378032 2. 00E -46 Unknown AK378032 
A109a BGIBMGA014427 0 Serine protease precursor AK378507 
A109b BGIBMGA000395 5.00E -77 Ribosomal protein L37 AK384780 


possibility is that the bacteria were removed from the 
body through feces (Shao et al., 2012). 

To identify the candidate genes that might be 
involved in the silkworm midgut response to infection, 
we employed ACP-based PCR technique to identify 
differentially expressed genes in the larvae after 
bacterial infection. Through sequence analysis, we 
identified 18 genes responsive to oral infection. This 
technique uses high annealing temperature during 
amplification and, thereby, significantly diminishes 
artificial positives ( Kim et al., 2004). Additionally , 
detection of PCR products on standard ethidium 
bromide-stained agarose gels reduces the cost. These 
benefits make the method a popular one for identifying 
differentially expressed genes in various studies. These 
include plant defense-associated genes ( Madala et al., 
2011) , and genes associated with ovarian cancer ( Kim 
et al., 2010) , metal metabolism (Han et al., 2010) , 
and host response to viral infection ( Balasubramaniam 
et al., 2011). This technique was also used in age 


estimation of forensically important blow fly ( Boehme 
et al., 2013 ) and insect immunity ( Freitak et al., 
2009; Bang et al., 2012). However, compared with 
microarray technique, this technique only makes a 
of differentially 
visualized and recovered on agarose gels. 

Insect gut immune system is composed of ROS 


small number expressed genes 


generation and antimicrobial peptides production. In 
the current study, we identified some genes that are 
known to participate in insect immune responses. We 
confirmed that B. mori PGRP-L1 responded to gram- 
negative bacterium PAOI infection. In Drosophila, 
PGRP-LC and PGRP-LE DAP-type 
peptidoglycan containing bacteria and induces the IMD 


senses the 


pathway ( Charroux and Royet, 2012). Manduca sexta 
PGRP-1 
initiates § _prophenoloxidase activation pathway 


(Sumathipala and Jiang, 2010). When infected by 
means of injection, B. mori PGRP-L1 found to be 


recognizes certain bacterial strains and 


induced by both gram-negative bacterium E. coli and 














































510 Be e+ Acta Entomologica Sinica 56 %& 
A B 16 
an P Pa 
E € 8 et Sa 
g z pag inan 
© ooh na Q AAE E ma 
O me $ 4 m pe 
E Bea is S = EES 
5 E 3 voce, Ep 
5) a tat oO u'n < ot 
= oe = =y f 
S T 5 l 
E em 
0.5 
1 2 8 16 24 1 2 8 16 24 
Hours post infection Hours post infection 
C D 
2 10 
Pa 
1 aa Sa 


0.5 


Fold (infection/control) 





0.25 


1 2 8 16 24 
Hours post infection 











aans 









0.001 


Fold (infection/control) 


0.0001 


1 2 8 16 24 
Hours post infection 


Fig. 3 Confirmation of four selected differentially expressed genes using quantitative real-time PCR 
A; Serine protease precursor gene; B; 30kP protease A precursor gene; C; PGRP-L1 gene; D; BGIBMGA004115 gene. Data in the figure are mean + SE. 


gram-positive bacterium S. aureus ( Tanaka et al., 
2008). Together with our finding, it appears that 
PGRP-L1 plays different roles in the midgut and the 
hemocoel. 

We also identified two serine proteases in the 
silkworm midgut, which responded to the infection. 
One of them ( BGIBMGA014427 ), a highly basic 
protease originally purified from the digestive juice, 
was thought to be involved in food digestion and 
cellular protein digestion during metamorphosis 
(Kotani et al., 1999). Our results showed that its 
gene expression increases after the silkworm was 
challenged with PAO1. 
(30kP protease A) gene was up-regulated remarkably 


The other serine protease 


upon infection. Serine proteases are key components of 
the prophenoloxidase activation pathway and Toll 


pathway ( Buchon et al., 2009 ). 


quantitative real-time PCR revealed that a group of 


Microarray and 


serine protease genes were up-regulated by bacterial, 
fungal and viral infection ( Zhao et al., 2010). 
Functions of these serine proteases in the silkworm 
immune system need further investigation. 

In summary, by using ACP-based differential 
expressed gene display technique, we identified 18 
genes that might be involved in gut immune response in 
B. mori. Expression patterns of four selected genes 
were confirmed by quantitative real-time PCR. PGRP- 
L1 and two protease precursor encoding genes were up- 
regulated by oral infection of PAOI and S. aureus. 
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